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PART I. NON-€QUILIBRIUM CALCULATIONS 
ABSTRACT 

The vibrational and rotational excitation of the CO molecule in cometary 
conae has been investigated using a model which includes IR vib*ational 
pumping by the solar flux, vibrational and rotational radiative decay, and 
collisional coupling among rotational states. Steady-state is not assuned in 
solving the rate equations. The evolution of a shell of CO gas is followed as 
it expands from the nucleus into the outer cona. Collisional effects are 
treated using a kinetic temperature profile derived from recent theoretical 
work on the ccma energy balance. The kinetic temperature is assumed to be 

extremely cold in the inner coma, and this has important consequences for the 
CO excitation. If optical depth effects are ignored, only low J transitions 
will be significantly excited in comets which are observed at high spatial 
resolution. Ground-eased observations of CO ro-vibrational and rotational 
transitions will be extremely difficult due to lack of sensitivity and/or 
terrestrial absorption, but CO should detectable from a large ccmet with 

favorable observing geometry, if the CO is a parent molecule present at the 
104 level (or greater) relative to Observations using cooled, 

space-borne instruments should be capable of detecting CO emission from even 


moderately bright comets. 



I. INTRODUCTION 


Although strong circvinstantial evidence Indicates that H^O ice is the 
dominant volatile constituent of ccmetary nuclei (Weaver ^ 1981, Delseime 

1982, Campins, Rieke, and Lebofsky 1983), ultraviolet (UV) observations of 
atonic carbon in cometary conae (summarized in Feldman 1982) clearly show that 
carbon- bearing compomds must com[»‘ise a significant fraction (>101) of the 
volatiles. CO is one likely source of this carbon. 

Observations have demonstrated that CO is present In some eometary comae 
although the CO abundance can not always account for the anount of atomic 

carbon observed. Feldman and Brtne (1976) observed several bands of the CO 

fourth positive group (A^-X^Z*^) in the UV spectrun of comet West (1976 VI). 

Feldman (1978) has analyzed these results and found that the CO production 
rate in this comet was •/‘301 of the H^O rate, and that the observed atomic 
carbon was derived from the CO. Fourth positive group emission was 
tentatively identified in International Ultraviolet Explorer ( lUE ) spectra of 
comet Brad field (1979 X) (Feldman et 1980 ). A*Hearn and Feldman (1980) 

demonstrated that the CO production rate (assuning CO to be a "parent" 
molecule) could be no more than of the H 2 O production rate in this comet. 
On the other hand, Weaver (1981) argued that atonic carbon was relatively 

abindant in comet Brad field (1979 X) (the production rate of the carbon parent 
was .rIOl of the H^O production rate), and could not be derived frcr-' either CO 
or COg. Cosmovici ^ (1982 ) recently claimed to have detected a triplet 

system of CO (d^A-a^n) in the visible spectrun of conet Brad field (1980 XV). 
Thus, CO is an important constituent of at least some cometary nuclei, even 
though it may not be the ultimate source of atomic carbcn in all comets. 
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Feldman (1983) haa suggested that the relative abisidance of CO may be one 
of the few fundanental characteristios which distinguish one ocmet from 

another. Such differences could provide clues to determine the formation 
regions of comets (Yanamoto, Nakagawa, and Fukui 1983). We have, therefore, 
Investigated the vibrational and rotational excitation of the CO molecule in 
canetary ccmae. The IR and rotational transitions of CO are potentially 
important probes of the density, kinetic temperature, and the outflow 
velocities of ccmetary parent molecules in the inner ccxna. 

Two recent papers have also addressed the problem of vibrational and 
rotational excitation of CO in comets. Krishna Swany (1983) has apparently 
neglected the most important vibrational excitation mechanism, direct 
radiative punping by absorption of solar IF radiation. Crovisier and Le 
Bourlot (1 983 ) have considered IR pumping but have concentrated on 
calculations in the outer cona where fluorescence equilibrium prevails, and 
they only briefly discussed the case of an isothermal inner coma where 

colllsional effects are important. They also pointed out that the time 
constants for the CO problem are long enough that steady state should not be 
assiwed at every point in the coma in solving the rate equations. 

In our approach to the CO excitation calculations, we have explicitly 
solved the time-dependent problem, i.e., we have followed the evolution of a 
shell of CO gas as it expands radially outward from the nucleus. In this 
manner, we have Included the inner cona where colllsional effects daninate the 

rotational excitation, the outer coma where radiative equilibrium dominates, 

3 5 

and also the intermediate region from 10 to 10 km from the nucleus where 
neither collisional nor radiative equilibrium is attained. In addition, we 

have considered a non-isothermal temperature profile for the inner coma. 

Details of the model calculations will be given in sec. II. From the model, we 
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have determined conditions under which the CO millimeter and IR transitions 
should be observable. Recently we have attmpted to observe sane of the 4.7 urn 
and Jrl«o CO transitions in comet IRAS-Araki-Alcock (1 983d). 

II. THE MODEL 


Before the model calculations can be performed, certain physical 
conditions in the comet must be defined. The kinetic temperature of the gas 
directly above the nucleus determines one of the initial conditions of the 
excitation model. Thermodynanlc considerations indicate that the temperature 

of a spherical, Isothermal nucleus composed principally of H^O ice is T^187K, 
for a comet at a heliocentric distance of 1 AU (Houpis and Mendis 1981). We 
adopt this value for the initial kinetic temperature of the gas leaving the 

nucleus. To calculate the effects of col'.isional excitation, the kinetic 
temperature profile of the inner coma must also be specified. Marconi and 
Mendis (1982) have demonstrated that, after release from the nucleus, the coma 
gas is rapidly cooled to extremely low temperatures (v5K at a distance of v^o 
km from the nucleus), followed by gradual photolytic heating beyond *r100 km. 
A figure in Marconi and Mendis (1 982) is used to derive a polynomial fit to 
their calculated kinetic temperature profile and this fit is the profile we 
use in our calculations. Gases subliming from the surface of the nucleus are 

assuned to flow radially outward with a velocity of 0.8 km s~^ (Weaver and 

Mumma 1984), Densities are calculated using a Haser model (Festou 1981), with 
lifetimes taken from Huebner and Carpenter (1979). H 2 O and CO are the only 
gaseous species considered, with H 2 O being dominant. As a result, collisional 
effects are confined to CO-H 2 O interactions. The H^O production rates are 
taken fV*om observations of comets, but these are generally consistent with 
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values derived from thermodynamic models. The CO production rate is taken to 
be iOt of the H^O production rate, but observed rates (or upper limits) are 
used for specific comets, whenever this is possible. We assume that the 
radiation field in the coma is determined solely by the solar field, with 
further details given below. 

Only the CO vsO and vs1 levels of the ground electronic state are 

considered in the present analysis (see discussion in Crovisler and Le Bourlot 
1983). Within £3ch vibrational level we use 21 rotational states since the 
results of the present work and of Crovlsier and Le Bourlot (1983) show no 
need for including rotational states higher than Js20. A Dunham expansion with 
coefficients from Kildal, Eng, and Ross (1974) is used to calculate the 
energies of the states, and the allowed transition frequencies are derived 
using these energies. 

The excitation of an expanding shell of CO gas is determined primarily by 
radiative pumping due to absorption of solar photons in r o-v ibrational 
transitions of the IR (1,0) band, radiative relaxation in the ro- vibrational 
lines of the IR (1,0) band and in the rotational lines of each vibrational 
level, and collisional coupling among the rotational states of each 
vibrational level. Collisional coupling of the vibrational levels is 
negligible throughout the coma (Weaver and Mumma 1984 ). It is easily 
demonstrated that stimulated emission in the ro>vibrational and rotational 
lines can be neglected, and that stimulated absorption in the rotational lines 
can also be emitted. The optical depths of seme low J lines (in both the IR 
(1,0) band and the rotational lines) may exceed unity close to the nucleus 
(r^lOOO km). However, considering the Increased complexity of a full 
radiative transfer treatment, we have assumed that all radiative transitions 
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•re optically thin. The problem of radiative transfer is, therefore, 
canpletely decoupled fron the rate equations. 

Thus, the rate equations for our model are given by: 


dCj/dt = 


1,0 1 1,0 1 0,0 0 

*J+1,j"j-»-1 ^ *J-1,j”j-1 


for vsO, and 


(1 ) 


dnj/dt 


- 0,1 0 „ 0,1 0 , 1,1 1 1 

J+1 ®J-1,J^"j-1 *^J+1,J"J+1 ^^k,A 

r * 1.0 »'»0 1 1 

■ ^^J,J+1 *'j,J-1 * ^J,J-1 ■*“ J 


for vsl. 


The rate equations given by (1) have the vibrational states indicated by 
superscripts and the rotational states InJ'cated by subscripts. When two 
superscripts or subscripts are used, the first refers to the initial state and 
the second to the final state. The vibrational index, v, can take values v=0 
and vsl, while the rotational index, J, runs from J sO to Js20. The unknown 
quantities are the nj*s which are the fractional populations of the rotational 
states of CO. 

Einstein A-coefficients for the ro-vibrational transitions (Nuth and Donn 
1981) are given by; 

P branch; s (64ir“/3h) (TM)^u^(J+1 )/(2J+1 ) 

s 3.39x10“^ij^(J+1)/(2J+1) (2) 

R branch: s 3.39x10”^u^ JA2J+1), 

where the transition frequencies (u) are given in cm'^ and the transition 
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monent (TM) is equal to 1.04x10 esu. Einstein B-coefficients are calculated 
in the usual way from the A-coefficients. We have expressed the radiative 
absorption rate constants in the IR (1,0) band as: 

P branch: = (F./8irc) (2J--i)/(2J+1 )a1’° , 

J,J-1 0 j-1,J 

s 33.2u'^(2J-1)/(2J+1 )a]^° J (3) 

R branch; s 33.2u"^(2J+3 )/(2J+1 )aJ^° j . 

The absorption rate constant in the IR (1,0) band is calculated assixning a 

13 - 2-1 -1 

constant solar flux over the (1,0) band of F^s2.5x10 photons cm s /cm at 
1 AU (Labs and Neckel "968). 

Rotational Einstein A's are given by (Somerville 1977): 

0 0 U R R '' 

A/, . = (64 7rV3hc^)v^u.'-J/(2J+1 ) 

J f J ••I W 

= 1. 164x10“^°v^Mq^J/(2J+1) (4) 

J-1 = 1.164x10 ^ J/(2J+1), 

with the transition frequencies (v) in MHz, and with dipole mcments (^^=0.1098 
and =0.0848 Debyes) derived from the constants given in the work of 
Kirschner ^ £l. (1977 ). 

The collision rate constants have the following form; 

-3 2 

where n^^^Q i’® ^2^ number density ( on ) , Oj j, is the cross section (cm ) 

for the transition from the J-th to the J '-th level of CO produced by 
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collisions with H 2 O, and v is the average relative speed (cm s*S of the H 2 O 
and CO molecules. 

The collision cross sections deserve some discussion because of the 

difficulties in estimating their magnitude. Crovisier and Le Bourlot (1983) 

point out that both experimental and theoretical data on this problem are 

-15 2 

scarce. In their calculations, they assune a cross section of o=2x10 cm 
(presumably this does not refer to a transition from one specific rotational 
state to another state, but to the total rate out of one rotational state), 
which they call the "geometrical" cross section. However, we calculate the 
geometrical cross section to be about an order of magnitude smaller than this 
(from Herzberg 1945, the separation C-0 in CO is 1.13 S while the separation 
H-0 in H 2 O is 0.958 fi) . 

A better estimate for the collision cross section than the geometrical 
value can be obtained by considering the data of Varghese and Hanson (1981), 
who have investigated the broadening of CO lines in the IR M,0) band by 
collisions witn H20* Using the theory outlined in Draegert and Williams 
(1968), we can use these data to derive a total cross section for the 
de-excitation and excitation of a CO rotational state by collisions with H 2 O: 

°tot ■ ■'•32x10“"'^ cm^. (6) 

Strictly speaking this cross section is valid only at room temperature 
(>t 294^98K), but these cross sections generally show little variation with 
temperature (S. Green, private communication). We have used this value 
throughout the coma. 

In addition, we must also know how collisions connect individual states. 
We assune that the cross section is a function of aJ only, with the explicit 
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dependence given by Table 1, This dependence on AJ is somewhat arbitrary, but 

was suggested to us by S. Green ([^ivate communication) , based on his 

experience with other molecules. By experimentation, we have found that our 

model is rather insensitive to the explicit aJ dependence of the cross 

sections, as long as the total cross section ronains fixed (compare this with 

results of H^-CO collisional excitation in Goldsmith 1972, and McKee ^ al , 

-14 2 

1982). We force the total cross section to always equal 1.32x10 cm . 
Excitation cross sections are related to de-excitation cross sections by the 
detailed balance condition (cf.. Goldsmith 1972). 

To complete the expression for the collisional rate constants, the 
relative thermal velocities ( v) are calculated assuming a Maxwellian 

distribution of velocities (Which is incorrect, but we neglect this fact) and 
the kinetic temperature profile of Marconi and Mendis (1982). When the 
collisional rates are larger than the radiative rates, then the detailed 

balance condition will ensure that the rotational excitation temperature is 

equal to the local kinetic temperature, i.e., rotational LTE will be satisfied. 

With all the necessary physical parameters of the comet model defined, 
equations (1) lead to 42 coupled, linear differential equations (LDEs) with 
only the collisional rate constants explicit functions of time. Since we 
assume that a shell of CO gas expands uniformly once released from the 
nucleus, we have. r=vt, where r is the distance from the nucleus, v is the 

expansion velocity (0.8 km s“^ ) , and t is the time. Therefore, a solution of 
our coupled LDEs will give the excitation of CO as a function of time, which 
is easily translated to distance from the nucleus. 

An IMSL ( International Mathematical and Statistical Libraries ) subroutine, 
DGEAR, ^.-iiich employs a Gear method (Gear 1971 ) for "stiff” equations, was used 
in sol/ing the system of coupled LDEs given in (1). A system is said to be 
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"stiff" when the time constants in the problem are of very different 

2 

magnitudes. Since the comet model involves rate constants ranging f« om vio 

“1 “8 “1 

s (collision rate constants near the nucleus) to *''10 s (rotational 

Einstein A for the J=1-*-0 transition), this numerical tecnnique is clearly 
necessp. y. Other techniques for solving LDEs, such as Hunge-Kutta methods, 
uer'i attempted on equations (1), but with little success. 

The initial conditions for equations (1) are obtained by setting dnj/dt=0 
at the surface of the nucleus, and solving the resultant set of simultaneous 
linear equations. The density and kinetic temperature are assimed to be 
values appropriate for r=2.5 km. Of the 42 equations of the homogeneous 
system, only 41 are independent, so we replace one of the equations by the 
normalization condition; 

}| = 1. (7) 

The 42 linear equations are then rioived by use of an IMSL routine LEQT2F, 

With the initial conditions as input, the DGEAR package attempts to 
integrate the coupled LDEs in time. As the program steps into different 
regions S the coma (or, equivalently, steps forward in time), rate constants 
are updated to take into account the density dilution due to expansion and 
lifetime effects, as well as the changing kinetic temperature. The DGEAR 
package autcnatically takes as small a step as necessary in order to give 
consistent results in adjacent steps. The solution consists of a set of 
relative populations for each energy level throughout the coma. We generate 
output on a logarithmic scale, using 10 points/ decad e (in km space) ranging 
from 10 km to 10^ km. The program was run on a Digital Equipment Corporation 
VAX 11/780, and a complete solution takes v'12 hourr of cpu time. 
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The relative fractional populations are used to calculate volune emission 
rates using: 


jJ*j,(r) 5 nj(r)Aj’j,n^Q(r)/(4ir), (8) 

V V* -'5 -1 -1 

where J,*.,(r) is the volune anission rate (photons cm ^ s sr ) for the 

V 9 J 

transition from initial state J to final state J’ at the position r, n (r) is 

wU 

the total number density of CO molecules at position r, and v and v’ are the 
initial and final vibrational levels respectively. Emission rates for the IR 
(1,0) band are calculated using v=1 and v’=0, while rotational transitions are 
calculated using vsv'=0. 

The surface (colunn) brightness is calculated by integrating j^’'[,(r) 
along a given line of sight. This integration is done numerically using 
path length that extends over 12 decades (in c:;. space; from 10° cm to 10 
cm), employing 18 points/decade. Since the line of sight integration requires 
volune emission rates for many more points than are output from our model, we 
use a logarithnic interp'lation procedure to fill in between the calculated 
points. For comparison with observations, we average the surface brightness 
over the field of view (FOV) of the instrument. 

III. RESULTS AND DISCUSSION 

Calculations have been performed using a variety of model input 
parameters. Our standard ’’bright" comet has an H2O production rate of 
® ® production rate of Q^q=2x 10 mol s” . Since the 

Marconi and Mendis kinetic temperature profile is strictly valid only at R 
(heliocentric distance) = 1 AU, we have concentrated on calculations for 
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comets at or near this value. However, in order to gauge the effects of 
fluorescent pumping, we have varied the solar flux to approximate the 
conditions of a comet at R=2 AU and R=0.5 AU (production rates and other 
parameters remain identical to their values at R = 1 AU). Using measured H^O 
production rates and upper limits on the CO production rate (Feldman, private 
communication), we '.‘.ave also considered the case of comet IRAS-Araki-Alcock 
(1983d), a comet which passed very close to the earth (minimum distance = 

0.032 AU) in the spring of 1983. For the week during which comet 1983d was 
observable from the northern hemisphere, its heliocentric distance was nearly 
constant with R=1 AU. Finally, we have considered three optimum viewing 

aspects during the apparition of comet Halley in 1935-1986. The complete set 
of model results can be found in Chin and Weaver (1984). 

Here, we concentrate on the most important effects which determine the 
excitation of CO in a cometary coma, and compare the non-equl ibriura 

s-alculations with those assuming steady-state at each point in the coma. The 
rotational excitation of the v=0 and v=1 levels are examined in detail as a 
function of distance from the nucleus. By neglecting the rotational structure 
we also elucidate several inportant points concerning the vibrational 
excitation. Finally, we average our results over different fields of view to 
determine which unique aspects of the non-equilibrium model would be 

observable. 

Figire 1 shows typical results from the non-equil ibrium model (hereafter 
referred to as the DGEAR model) for the standard bright conet. A single 
contour in Figure 1 represents a fixed relative population. The contour 
levels are given in logarithmic units so that, for example, the curve labeled 
-1.0 represents a relative population of 10%. The vertical scale is the 
distance from the nucleus, while the solid and dotted vertical lines represent 


in 

the J values for the various energy levels. The regions between the vertical 
lines have no physical significance. Similarly, Figure 2 gives the contour 
plot for the excitation temperature, T , calculated between adjacent 

OX 

rotational states within a vibrational level 

(T =(Ey ,-Ey)/k/lnC(ny/ny ,) (2J-1)/(2J+1 )]) , vrtiere E is the energy of the 
state. Relative populations and excitation temperatures are plotted for both 
the v=0 and v=1 levels. The results are shown in this way because the contour 
plot is a convenient way to present data obtained for all states 
simultaneously. 

Several important points are apparent upon examination of the contour 

plots. In general, the results of the calculations behave in the expected 

manner. For r<10 km, the rotational states are maintained in an extremely 

cold distribution; i.e., only the lower J states are substantially populated. 

In this region, the H^O density is high and the rotational excitation 

temperature is expected to follow, at least approximately, the local Kinetic 

temperature frofile (however, see discussion below). In the outer coma (r>^10^ 

km) the rotational distribution is spread over many states as the molecules 

approach fluorescence equilibrium. For practical purposes, a fluorescence 
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equil ibri'-4n distribution can be assumed whenever r>10 km. The region between 

3 0 

10 and 10 km from the nucleus is a transition region in which radiative 
effects begin to take over from collisional effects. 

Fig ires 3 and 4 illustrate more clearly the behavior we have outlined. 
These show the fractional populations for the JsC rotational state of the v=0 
and v=1 vibrational levels, respectively. The results for the DGEAR model 
output, a single cut along the J=0 line in Figure 1, and for a homogeneous 
solution, which assumes steady-state throughout the coma, are plotted together 
in Figures 3 and 4, The solid curve represents the fractional population if 
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local theraodynanic equilibrium (LIE) is achieved ev«*ywhere. In general, the 
fractional population of the vsO level follovs an LTE distribution in the 
inner coma as can be seen in Figure 3. Both the DGEAR and homogeneous model 
results closely follow the solid curve for r<lO^ km. (towever, as Figure 4 
shows, rotational LTE is never achieved anywhere in the coma for the JsO, v=1 
state. In fact, rotational LTE is not achieved anywhere in the coma for any 
rotational state in the vri level. We will discuss this point in more detail 
when we examine vibrational excitation later. 

We also see from Figure 3 that for the v=0 level the DGEAR result departs 

from LTE at a much smaller radius (r^50 km) than does the homogeneous 

solution (rjJiOOO km). This is to be expected for several reasons. First, the 

HgO density is decreasing at increasing distances from the nucleus, implying 

that the collisional rate constants, which determine the equilibration time, 

are also decreasing. In addition, the kinetic temperature profile changes 

dranatically in this region of the ccma, reaching a minimixn at *^0 km but 

heating up at larger radii. Since the collisional rate constants are 

- - 1/2 

proportional to the square root of the temperature (C, ,,“V and v«*T ), the 

rate constants will also exhibit a fairly large variation. This combination 

of small and rapidly changing rate contants means that the DGEAR solution 

cannot adjust to the changing kinetic temperature and, thus, lags behind the 

homogeneous solution, which assixnes that steady-state has been achieved. In 

effect, the DGEAR solution retains a "memory” of the coldest portion of the 

coma, even as the coma tanperature increases. The populations in the DGEAR 

solution change in a smoother fashion than the homogeneous result, especially 

3 4 

in the intermediate region (10 -10 km) when the homogeneous solution 
continues to adhere to LTE, 


A further consequence of the smooth departure of the DGEAR solution from 
an LTE distribution is illustrated in Figure S, which shows the excitation 
temperature for the Jsl-K) transition in the vsO vibrational level. Again we 
plot a DGEAR solution, a homogeneous solution, and a solid curve which 
refx'esents the assumed kinetic temperature profile. The homogeneous solution 
shows an excitation temperature which remains thennalized to km, as is 

expected by comparing the collisional rate constants with the rotational 
Einstein A-coefficients (the disequilibrium distance in Weaver and Mumma 
1984 ). However, the DGEAR solution departs from the kinetic temperature at a 
much smaller distance from the nucleus C r^'TOO km) and becomes highly 
suprathermal , even inverted, in the transition region. Once again, this is a 
consequence of the '’relaxation” required by the DGEAR model as the 
distribution evolves from LTE to fluorescence equilibrium. While the J=1-*0 
transition becomes suprathermal or inverted, all other transitions depart from 
LTE by becoming subthermal for the DGEAR solution. 

A DGEAR solution using an isothermal (T=20CK) temperature profile with the 
physical parameters of the standard bright comet does not show any temperature 
inversions in the transition region. This result is expected since the 
fractional populations of a 200K and fluorescence equilibrium distribution do 
not differ greatly; essentially no "relaxation” is necessary to connect the 
thermal ized inner coma to the radiatively determined outer coma. It is the 
combination of a rapidly changing temperature profile and the use of the DGEAR 
method v*iich results in the large deviations from thermal ized populations and 
excitation tanperatures . 

Temperature inversions occur for several rotational states in both the v=0 
and vs1 levels when other model parameters are used (see results in Chin and 
Weaver 1984). Since the rotational Einstein A's increase as J^, there is a 
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tendency to form a bottleneck in the lower J levels, which have low decay 
rates (especially in CO which has a snail dipole mcnent). Inversion will 
result if no competitive excitation mechanism, such as collisional excitation 
or radiative pumping, is strong enough to distribute the rotational population 
among a nunber of J levels. For example, in the case of comets with low gas 
production rates or at large heliocentric distances, the lower J states in CO 
can easily become suprathermal or inverted. 

No appreciable anplification is expected from the inverted regions in the 
coma. An upper limit to the expected gain is <0.01% in our bright comet 
model. Indeed, we can continue to assune that the coma is optically thin. The 
excitation problem remains decoupled from the radiative transfer problem and 
the emergent radiation from the cona can be calculated in the usual manner 
(i.e., use equation 8). 

The relative populations of the upper and lower (v=1 and v=0 ) vibrational 
levels are determined only by the solar flux (collisional excitation is 
neglible between vibrational levels, as discussed earlier). If we ignore the 
rotational structure for the moment, statistical equilibrium gives the 
following steady-state vibrational distribution; 

nVn° = Bfl/A, (9) 

where A and B represent the band values for the Einstein A and B coefficients 
and p is the density of the solar radiation field. This solution represents a 
system vAiich is in equilibrium with the solar radiation field and fixes the 
ratio, nVr°, at 7.28x10”^ when Fs1 AU. The sun at a distance of 1 AU has an 
"equivalent” blackbody temperature of 261K at «r4.7 um, and that is the 
vibrational excitation temperature given by equation (9). The output of our 
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mod«l is consistent with this value and forms an important check on the 
correctness of the calculations. 

The explicit form for vibrational excitation as a function of time is 
given by (see the Appendix): 

nVn^ s Bp/A ♦ (10) 

where is a constant determined by the initial conditions. We see that the 

steady- state vibrational distribution is achieved very rapidly. The time 

constant for equilibration is 1/A (or 0,029 s) , not 1/Bp (or M.lUxlO s) , as 

might be expected. For t>>1/A, equation (10) gives the steady-state value. 

It is illustrative to exanine the case when only the lower state is 

initially populated (nVn^=0 at t=0 ) and for t<<1 /A. Then equation (10) 
1 0 

reduces to n'/n =Bpt, and the vibrational excitation is determined initially 
only by the solar pump rate constant (Bp). However, once some population is 
achieved in the upper level, radiative decay quickly determines the 
equilibrium value. This is due to the fact that since the radiative decay 
rate constant is so large canpared to the punp rate constant, the equilibrium 
value must be very small (i.e., Bp/A«1). 

Since the relative vibrational populations reach equilibrium almost 
instantly after release from the nucleus, the band excitation factor, or 
^band* fixed at a value of 2.^x10 s throughout the entire coma for Rsi 
AU. 

The intensities of individual transitions within the IR (1,0) band are 
determined by the rotational distribution of the v=1 level. As discussed 
earlier, the vsl level is never in rotational LTE anywhere in the coma. 
Figure U shows the fractional population of the J=0, vsl state for the DGEAR 
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solution, the hcmogeneous solution, and the LTE distribution. Both the DGEAR 
and homogeneous solutions give similiar results in the inner coma, an 

essentially flat distribution which is depressed by a factor of 2.5 - 3 from 
the peak of the LTE distribution at r</'100 km. Figure 6 shows the excitation 
temperature for the Jsl-^ vs1 transition for the DGEAR and homogeneous models, 
as well as the assumed kinetic temperature profile. Departure from the 

kinetic temperature profile is most pronounced at r>/'100 km, and both the DGEAR 
and homogeneous solutions give similiar suprathermal results throughout the 
inner cona. 

By examining equations (1) separately for the v=1 and vsO cases, we can 
understand why LTE is maintained in the v=0 level at the same time LTE is 

violated in the v=1 level. It is important to recognize the difference 

between the rate constants and the rates in equations (1). The rate constants 
are the Bp's, A's, and j,'s which multiply the fractional populations. The 
rates are the products of the rate constants and the fractional populations. 
In order to determine the dominant excitation processes for each vibrational 
level, it is necessary to detrnnine the dominant rates in (1). 

Although the collisional rate constants can be rather large (''50 s”^ at 
rslO km) with respect to B^'^p s“b and A^*^ (j<10-33 s“b in the inner 

coma, for the v=1 states the collisional rates (zCj j,n^) are much smaller 
s“b and must be compared to B^'^pn*^ s“^) and A^'^n^ 

s”S . We see, therefore, that qualitatively the collisional and radiative 
rates are about equal in the v=1 level. When the collisional rates decrease, 
for example when the kinetic temperature or the HjO density decreases, the 
radiative rates dominate. The largest departures from rotational LTE in the 
vs1 level occur in the coldest regions of the inner coma. On the other hand, 
in the vsO level rotational LTE is maintained within the inner coma. 
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The usual formula for the g- factor can be derived from equations (1) by 
assuming that radiative rates donlnate the excitation. If we neglect all 
colllslonal and rotational radiative decay terms, and assume steady-state, we 
get: 


[B ' pn. +B ’ ,pn. .]/(A * +A/, 

J+1,J J+1 J-1J J-1 J,J+1 J,J-1 


( 11 ) 


The g- factor is obtained simply by multiplying (11) by the .opropriate 
radiative decay rate. For example, is given by (with the transitions 

properly identified): 


®R0 ■ 


(Bp)p2n2]tARo/(ApQ+Ap2)]. 


( 12 ) 


The g-factor can then be evaluated by assiining a rotational distribution for 
the vsO level. Comparing the values of g^^ obtained by evaluating (12) 
assuning a rotational LTE distribution for v=0, to those obtained directly 
from the results of the DGEAR model, we get g .(DGEAR )/g (radiative)sl . 11 at 

nO n 0 

r = 10 km. Since g^^ (DGEAR)/gpQ(radiative) is also equal to 

n| (DGEAR)/n^^ (radiative) , we find that the radiative formula, equation (11), 
gives very good approximations to the actual fractional populations even at 
rslO km, where collisional effects are strong. 

The rotational states in the vzl level retain some information about the 
rotational distribution of the vsO level through equation (11), Although it 
appears that the vri level is in rotational LTE within the inner coma (the 
general trend of the contour lines for both v levels in Figure 1 are quite 
similiar), in fact, the upper v states are merely reflecting the LTE 
distribution of the vsO level. 
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Many of the differences between the XEAR and homogeneous models are 
avers^ed out by the line of sight Integration used in calculating column 
brightnesses. For example, Figure 7 shows the P1 colimn brightness versus the 
projected distance from the nucleus for both the DGEAR and homogeneous models, 
using the standard bright comet parameters at Rsl aU. The CO column density 
Is also plotted on the same graph. The colimn density falls off inversely 

5 

with the impact parameter until ^10 km when photolysis becomes important. 
The PI colimn brightnesses in Figure 7 differ significantly only when the 
impact parameter is j'10‘^-10 km (i.e., when the line of sight is about tangent 
to the transition regions in the cona) . The maximum differences between the 
two models are »WiOi. For higher J transitions, differences between the DGEAR 
and homogeneous results are also found in the inner cona, in addition to the 
transition region, with maximum differences reaching ^/50{ (these results can 

c: 

be found in Chin and Weaver 1984). At a distance of j' 10 km, where 

fluorescence equilibrium is achieved, colunn brightnesses for all transitions 
become proportional to the column density, as expected. 

Figure 8 shows colimn brightnesses for the PI and P7 transitions using 
only the DGEAR model. The P7 brightness is essentially flat for impact 

pcr'^ameters <10 km with a slight increase close to the nucleus. This behavior 
is indicative of the cold rotational distribution in the inner coma where 
higher J states are not significantly populated and, therefore, contributions 
to the line of sight integration can only come from outside of this region. 
The reverse is true for the PI transition for which the inner coma contributes 
significantly to the intensity. 

The intensity distribution in the IR (1,0) band changes dramatically for 

different fields of view. Table 2 gives aperture- aver aged absolute 

-•2 —1 

intensities (in W cm"* sr" ) for lines in the IR (1,0) band for two fields of 
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view (1* and 4") v^en our standard bright comet is 1 AU from the sun and 0.62 
AU from the earth. The detectability of these lines will be discussed in the 
next section. Flgtre 9 shows the aperture- aver aged g- factors (Which are 

A 

directly proportional to the intensities) for the 1* FOV. At a geocentric 

^ 4 

distance (a) of 0.62 AUt a 1' aperture subtends a radius of 1.35x10 km. We 

plot the results for the DGEAR and hcmogeneous models, and, for comparison, an 

LTE distribution using the assumed kinetic temperature profile shown in 

Figures 4 and 6. For this FOV and observing aspect, a large portion of the 

coma Is averaged into the measurement, and the intensity is distributed among 

many J states. For the DGEAR and hcmogeneous models, the peak intensity is 

reached at J=6 or 7, while the LTE distribution is peaked at even higher J and 

exhibits a "flatter" spectral profile. The "hot" LTE distribution is due to 

the assumption of high kinetic temperatures (T*/*500K when r>10 km) in the 

J' 

outer cona. 

In contrast, Figure 10 shows the g-factors for the IR (1,0) band averaged 

A 

over a 4" diameter circular aperture at the sane heliocentric and geocentric 
distances. Here the projected radius is 900 km, so that the measurement is 
primarily sampling the inner ccma. In this case, the DGEAR, hcmogeneous, and 

A 

LTE distributions all give similiar results. Note that for a 4" aperture, the 

A 

peak g-factors are almost three times as strong as the peak values for a 1' 
aperture. Only the low J transitions are significantly excited in the cold 
inner ccma, so that almost all the intensity in the band is concentrated into 
these few lines. The differences between the model calculations and the LTE 
distribution result from the fact that the v=1 rotational states are not in 
LTE (as discussed above) . 

Previous studies have obtained similiar values for the band g- factor 

-4 -1 

band"^'^ ® ’ while Weaver and 


(Crovisier and Le Bourlot 1983 derive g. 
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Munma 1984 give 2.5x10 s , canpared to 2,4x10 s frcxn our model 
calculations). The small differences arise from the different values used for 
the Einstein B-coefficients and for the solar flux. However t the earlier 
works make very different assumptions concerning rotational excitation in the 
coma, and, thus, derive different intensity distributions within the IR (1,0) 
band. Crovisier and Le Bourlot concentrated on a fluorescence equilibrium 
distribution and their results are simillar to ours for comets which are 
observed with large apertures at large geocentric distances. Weaver and Mumma 
asstned an LTE distribution for the v=0 level with T=200K; their results are 
similiar to those obtained assuming fluorescence equilibrium. For small 
fields of view, our calculations become critically dependent upon the assuned 
kinetic temperature profile in the inner coma. The "cold" IR intensity 

A 

distribution predicted for observations with a 4" FOV is a reflection of the 
extremely cold kinetic temperature profile used in our calculations. 

Table 3 shows the intensities of the rotational transitions for the bright 

A A 

comet model at R=1 AU a=0.62 AU for both 1' and 4" fields of view. The 
intensities in Table 3 are calculated assiming a line width of 1.6 km s”^ and 
are expressed as brightness temperatures in the Rayleigh-Jeans limit. The 
Crovisier and Le Bourlot (1983) results are not directly comparable to ours, 
since there are mmerous differences between the two models. However, if we 

A 

assure the same conditions as Crovisier and Le Bourlo- (i.e., 0.5’ FOV, A=1 
AU, and use the same Haser model input parameters and cometary line width) , we 
calculate significantly different intensities for some of the rotational lines 
(the largest difference is 1.8 for the Jsl*0 transition). Again, the higher 
spatial resolution results in Table 3 show enhanced values for the low J 
transitions and are a reflection of the "cold" kinetic temperature profile we 
have assuned for the inner coma. 
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IV. OBSERVATIONS 

The preceding discussion has important implications for the detectability 
of canetary anission in the IR (1,0) band and from the millimeter-wavelength 

A 

rotational transitions. The strongest IR line in a 4” FOV when ^=0.62 AU is 

expected to be the P2 transition with an average brightness of 1.6x10 Watts 

-2 -1 

cm sr (see Table 2). The sensitivity of the Kitt Peak National 

Observatory (KPNO) Fourier Transform Spectrometer (FTS) at the 4-meter Mayall 

Telescope in the region of the CO IR (1,0) band, for a one hour integration 

-1 -11 -2 

period, using its highest resolution (0.0216 cm ), is >^x10 Watts cm 

sr“^ (1 a limit) . This sensitivity limit can be achieved for cometary 

observations when they are performed at low air mass (^1), and if the cometary 
line is Do ppl er- shi fted far enough away from the peak of telluric CO 

absorption that the transmittance at the position of the cometary emission is 
essentially the continuum transmittance. Cometary geocentric velocities 
greater than 30 km s”^ are necessary in order to accomplish this, but such 
velocities are not unusual. Thus, it appears possible, although difficult, to 
do ground-based observations of the CO IR (1,0) band emission from a bright 
comet if the CO abundance is >10% of the H^O abundance. 

Ground-based FTS observations are background noise limited in this 
spectral region. Even when a narrow filter is used to limit the bandwidth to 
v'100 cm”\ covering most of the CO (1,0) band, background noise is still about 
two orders of magnitude above detector noise limited sensitivity. On the 
other hand, an FTS on a cooled, spac^borne instrument, like the Shuttle 
Infrared Telescope Faciltiy (SIRTF) , should easily detect cometary CO emission 
in the (1,0) band. 


Va 


25 


Ground-based searches for rotational transitions will be confined 

primarily to low J transitions like and since higher J 

transitions are difficult to observe from the ground and no sensitive 

receivers are currently available for transitions (see Wilson 1983). 

Since tjplcal millimeter-wavelength telescopes have rather large beam sizes, 
the relevant intensities to examine in Table 3 are the average values in a 1 ' 
FOV. Millimeter-wave receivers can now achieve sensitivity limits approaching 
v/'IO nK do) for the J=1-^0 transition. As can be seen from Table 3, the 

highest possible sensitivities are required to observe this transition from 
the ground for a bright conet. Receiver sensitivities are not as good for the 

J=2-f1 transition but since this transition will generally be brighter in the 

comet, it would not be unreasonable to also search for this line. The 

prospects of detecting either of these lines drop dramatically when the 

comet's geocentric distance becomes large (A>1 AU ; see earlier discussion). 
Since the CO brightness peaks strongly towards the nucleus if CO is a parent 
molecule, accurate ephemerides are essential for centering the bean on the 
nucleus to improve the sensitivity of any observation. 

The apparition of comet IRAS-Araki-Alcock (1983d) in May 1983 provided a 

favorable opportunity to perform sensitive IR and radio searches for CO. We 

searched for individual lines of the IR (1,0) band on 198? May 10 

( »r1 0: 20-1 1 : 31 UT) using the KPNO FTS at the Mayall 4 -meter telescojpe (These 

observations are sunmarized in Weaver ^ al. 1983). We estimate that the 

atmospheric transmittance for the strongest CO line was >^0.25 (the comet was 

observed at v«2 air masses). Ito CO lines were detected. Our upper limit do) 

1 4 «2 

for the average CO colunn density in the 4” FOV is 1.3x10 cm , which 

27 -1 

translates into a CO production rate of v'lxlO mol s . From IDE 


26 


observations of OH anlssion, the H 2 O production rate is estimated to be 
between 1.5-4x10^® mol s“\ 

A search was conducted for the CO Js1-^o radio transition on 3 nights C983 

May 7f10» and 11 UT) using the Bell Laboratories 7-metei telescope in Holrodel, 

New Jersey (These observations are also described briefly in Weaver ^ 

1983). Telescope poin<-^ng and tracking were verified using a TV guider system 

that was co-aligned with the radio telescope. The small geocentric distances 

(a^O. 07-0.03 AU) of the comet ensured that the coldest region of the cona was 

observed (the beam size was .rlOO ") , and that there would be reasonable 

expectations for observing the J=1-»0 line. Unfortunately, no CO emission was 

detected. The upper limit on the CO production rate based on the radio 

results is virtually identical to the number derived above from the IR 

27 -1 

observations (e.g., Qj,q<1x 10 mol s ). 

This upper limit is also consistent with the upper limit placed on CO 

production from lUE observations (fourth positive group emissions were not 

P 1 

seen in the IU£ spectra). The lUE limit (Q^q<4x 10 mol s” on 5/12 UT; 
Feldman, private communication) is more sensitive than the IR and radio limits 
but, unlike the IR and radio observations, the lUE measurements would yield no 
direct information on the rotational excitation or outflow velocities of the 
CO molecules (the lUE resolution is not high enough to resolve individual 
lines within the CO bands). If CO were a parent molecule present at the same 
abundance relative to H^O as for comet West (1 976 VI) (Q_yQ„^j«0. 3) , then CO 
emission should have been detectable in all 3 spectral regions. 

V. CONCLUSION 


Our study takes two previously-neglected factors into account in 
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calculating the vibrational and rotational excitation of CO. First » we 
explicitly derive the fractional populations of CO as solutions to 
time- dependent differential equations, so that the physical parameters of the 
comet determine the evolution of the solution throughout the coma. Important 
differences between the DGEAR method and a steady-state solution occw , 
especially in the transition region (rvIO to 10 km). The DGEAR solution 
predicts suprathermal and inverted populations for low J rotational 
transitions. Second, we consider a non- isothermal and cold kinetic 
temperature profile for the inner coma (r<10-* km). This leads to enhanced 
emissions for low J transitions in the IR (1,0) band and for 
millimeter-wavelength rotational transitions in the inner coma. 

The predictions in this paper are crucially dependent upon our choice for 
the kinetic temperature profile. If the kinetic temperature is very different 
from that given by present models (e.g., Marconi and Mendis 1982, Shimizu 
1976, Crovisier 1 983 ), then many of our conclusions concerning the 
observations of comets using snail fields of view or at small geocentric 
distances may be invalid. We have also neglected optical depth effects. 
Significant optical trapping could affect the excitation of CO in the inner 
coma, and this could alter our results. To our knowledge, a proper treatment 
of multi-level radiative transfer in comets has never been attanpted, and is 
outside the scope of the present paper. Such a refined analysis of the CO 
problem may not be justified until suitable future observations of CO indicate 
that this approach is warranted. 

The predicted intensities for our bright comet model show that IR (1,0) 
bind and millimeter-wavelength observations of CO are feasible under certain 
favorable conditions and would be important probes of the excitation 
conditions in the comae of comets. The kinetic temperature of the inner comae 


of ccxneti is expected to be cold and IR and mill line ter-wav el ength observations 
of CO would directly test this hypothesis. Due to msxiy factors, including 
terrestrial absorption, bean dilution, and lack of sensitivity, ground-based 
observations may only be accomplished for large comets passing reasonably 
close to the earth. Conet IRAS-Araki-Alcock (1983d) (u-esented an excellent 
opportunity to test many of the predictions of our model calculations, since 
it passed within 0.032 AU of the earth. Unfortunately, 1983d proved to be a 
snail comet with relatively low gas production rates (Large comets could have 
H^O production rates a factor of 10 or more larger than comet 1983d at the 
same heliocentric distance). 

Our calculations also indicate that comet Kalley will not be a favorable 

target for ground-based observations of CO. During the second week of March 

in 1986 (one of the most favorable times for observing comet Halley; R=0.88 AU 

and 1 AU), parent molecular colunn densities may be lower by a factor of <^7 

29 —1 

compared to those of comet 1983d (We assume ^ comet 

Halley when R=0.88 AU), Moreover, the rather large geocentric distance will 
mean that ground-based observers will be viewing a CC population closely 
approachi;i6 fluorescence equilibrium, and this will further inhibit the 
detection of the low J IR (1,0) band and millimeter-wavelength transitions. 

Although ground-based IR observations of CO missions from comets will be 
extremely difficult, a completely different picture emerges when one considers 
observations fron earth orbit using cooled instrunentation. An FTS instrument 
aboard a cooled telescope such as SIRTF could have up to two orders of 
magnitude increased sensitivity canpared to similiar ground-based instriments 
in the CC IR (1,0) band. Thus, Earth orbiting instruments should be capable 
cf monitoring CO emission frcm even moderately-bright comets. 
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APPENDIX 

We wish to solve the explicit time development for a 
which is excited only by radiative processes. Let the upper 
by n^ and the lower state by n^. We define R as 

R 3 n,/i^. 

Fran the normalization condition, n^+n^sl, we have: 
dn^/dt = -dn^/dt. 

The time derivative of R is given by: 


dR/dt - (3R/3n^)(dn^/dt) + (sR/sn^KdCQ/dt) 

= (1/n^)(dn^/dt) - (n^/np)(dnQ/dt) 

- -(1 /nQ)(dnp/dt) - (R/njjXdn^/dt) 

s -(dnQ/dt)(l+R)/nQ, 

using equations (A1) and (A2). The rate equation is: 


dCg/dt = An^ 


- Bpn, 


• ‘""o • “""o 

= n^CAF - Bp) , 


two level system 
state be denoted 

(A1) 


(A2) 


(A3) 


(A4) 
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where A Is the radiative decay rate constant and Bp is the radiative pimp rate 
constant. StA>stitution from equation (A4) into (A3) yields: 

dH/dt = Bp ♦ (Bp - A)R - AR^. (A5) 

The equilibriun solution can be obtained by setting dR/dtsO. and solving the 

resultant quadratic equation. The roots of equation (A5) are RsBp/A, and Rs-1 

(this gives a nor>-physical solution). Setting equation (A4) to zero and 

solving for R also gives the equilibrium result RsBp/A. Dividing equation 

2 

(A5) by -A, and ignoring terms in R and (Bp/A)R (R<<1 and Bp/A<<1) gives: 
-(l/A)dR/dT = R - Bp/A. (A6) 

Equation (A6) can be rewritten as 

dR/(R-0p/A) = -Adt. (A7) 

Equation (A7) can be Integrated inriediately to yield: 

ln(R-Bp/A) = -At + C, (A8) 

which can be exponentiated to give; 

R r Bp/A + R^e"^^, (A9) 

where R^ is determined by the initial conditions. Note that the time-constant 
in equation (A9) is dependent only on A, and not on the pumping rate. 
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TABLE 1 

CO-HgO Collision Cross Section Information^ 


aJ (=J -J, ) Fraction of Total De-Excitation 

upper lower 


1 

2 

3 

4 

5 

6 
>6 


0.34 

0.25 

0.20 

0.10 

0.07 

0.05 

0 


a -14-2 

The total cross section is always forced to equal 1.32x10 cm . 

Excitation cross sections are derived from de-excitation cross sections by 


using the detailed balance relationship 
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TABLE 2 


Predicted Intensities for the CO IR (1,0) Band; Bright Cornet^ 


Transition 


Frequency (cm 


) 


Intensity (10**^^ W 


-2 - 1 . 
cm sr ) 


1 • FOV 


4*' FOV 


P12 

209‘«.863 

0.162 

0.344 

P11 

2099. 083 

0.209 

0.503 

P10 

2103.270 

0.256 

0.732 

P9 

2107. '124 

0.298 

1.08 

P8 

2111.543 

0.333 

1.62 

P7 

2115.629 

0.355 

2.51 

P6 

2119.681 

0.365 

4.00 

P5 

2123.699 

0.361 

6.43 

P4 

2127.683 

0.343 

10.3 

P3 

2131.632 

0.310 

15. 1 

P2 

2135.547 

0.241 

15.8 

PI 

2139.427 

0. 131 

9.41 

RO 

2147.082 

0.123 

8.06 

R1 

2150 . 856 

0.214 

10. 4 

R2 

2154.596 

0.270 

8.12 

R3 

2158. 300 

0.308 

5.49 

R4 

2161.969 

0.329 

3.61 

R5 

2165.601 

0.335 

2.36 

R6 

2169.198 

0.324 

1.58 

R7 

2172.759 

0.300 

1.08 

R8 

2176.284 

0.264 

0.755 

R9 

2179.772 

0.221 

0.532 

RIO 

2183.224 

0.175 

0.372 

1 production 

P R 

rate is 2x1 0 

, -1 
moi s 

and the H-0 production 


is 


29 

2x10 mol 


-1 

s . 


The heliocentric distance is 1 AU and the geocentric distance 


is 0.62 AU. The average CO colunn density in the 1’ circular field of view 
(FOV) is 8.55x10^^ cm"^, while that in the 4” circular FOV is 1.38x10^^ cra"^. 
All intensities are aperture- aver aged values calculated using the 
non-equilibrium (DGEAR) model. 
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TABLE 3 

Predicted Intensities for the CO Rotational Transitions: Bright Comet® 


Brightness Tanperature*^ (K) 


Transition F requency (GHz) 1 * PQV 


4” FOV 


1-0 

115 

0.008 

0.595 

2-1 

231 

0.027 

1.18 

3-2 

346 

0.050 

1.28 

^♦-3 

461 

0.076 

1.15 

5-4 

576 

0. 101 

0.940 

6-5 

691 

0.123 

0.746 

7-6 

807 

0. 138 

0.588 

8-7 

922 

0.144 

0.465 

9-8 

1037 

0. 141 

0.369 

10*9 

1152 

0.128 

0.289 


3 ^ 8 “ 1 

The CO production rate is 2x10“ mol s” and the H 2 O production rate is 
29 -1 

2x10 mol s . The heliocentric distance is 1 AU and the geocentric distance 
is 0.62 AU. The average CO colunn density in the 1' circular field of view 

(FOV) is 8.55x10^^ cm"^, while that in the 4*' circular FOV is 1.38x10^^ cm”^. 
All intensities are aperture- aver aged values calculated using the 
non-equilibrium (DGEAR) model. 

*^The Brightness Temperature is calculated in the Rayleigh-Jeans limit assuming 
a rectangular line width of 1.6 Icu s 
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FIGURE CAPTIONS 

Figure 1. Contour plots showing relative populations of CO rotational states 
as a function of distance from the nucleus. Each vertical line represents a 
different rotational state, starting with J =0 at the extreme left and 
incrementing by one until J=20 is reached at the extreme right. The regions 
between the vertical lines have no physical significance. The contours arc 
given in logarithmic units so that, for example, the curve labeled -1.0 
represents a relative population of 10J. Plots are shown for both the v=0 and 
vs1 vibrational levels of the ground electronic state. The non-equil ibrium 
(DGEAR) excitation model is used to generate these results. The CO production 
rate is 2x10 mol s , the H^O production rate is 2x10 ^ mol s” , and the 
heliocentric distance is 1 AU. 

Figure 2, Contour plots showing the excitation temperatures for CO rotational 

transitions as a function of distance from the nucleus. Excitation 

temperatures are calculated between adjacent rotational levels (see text for 

explicit formula) with each transition represented by a vertical line. 

Excitation temperatures for the lowest frequency transition (J=1*0) are given 

along the left most dotted vertical line. Notice that this transition for the 

v=0 level becomes inverted (negative excitation temperatures) in the region 

between 3000 and I 300 O km from the nucleus (this region is cross-hatched in 

the figure) . Plots are shown for both the v=0 and v=1 vibrational levels of 

the ground electronic state. The non-equilibrium (DGEAR) model is used to 

28 —1 

generate these results. The CO production rate is 2x10 mol s” , the H^O 

29 -1 

production rate is 2x10 mol s , and the heliocentric distance is 1 AU. 
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Figure 3. The frcctlonal population of the JsO, v=0 level of CO is plotted as 
a function of distance from the nucleus. Three different profiles are shown 
corresponding to thret different excitation models: the circles, 0, represent 
results from the non-equlibrium (DGEAR) solution, the triangles. A, show the 
steady-state results, and the solid line is the profile if rotational LTE is 
valid eve’*ywhere in the coma. The CO production rate is 2x10 mol s” , the 
HjO production rate is 2x10“ mol s , and the heliocentric distance is 1 AU. 

Figure 4. Same as Figure 3, except for the J=0, v=l level of CO. 

Figure 5. The excitation temperature of the J=1+0, v=0 transition in CO is 

plotted as a function of distance from the nucleus. The circles, 0, show the 
non-equilibrium (DGEAR) solution, the triangles. A, are the steady-state 
results, and the solid curve is the assumed kinetic temperature profile. For 
the DGEAR results there is a region in the coma (3000 km < r < 13000 km) where 
the excitation temperature becomes negative (indicating a population 
inversion). For these points we plot the absolute value of the excitation 

temperature, and indicate that ’nversion has ocurred by arrows. The CO 

38 "1 29 “1 

production rate is 2x10 mol s , the H2O production rate is 2x10 ^ mol s” , 

and the heliocentric distance is 1 AU. 

Figure 6. Same as for Figure 5, except for the J=1-0, v=1 transition in CO. 
In this case, there are no population inversions, but both the non-equilibrium 
(DGEAR) and steady-state results indicate suprathermcl excitation temperatures 
in the inner ccma. 


Figure 7. Colunn densities and colttrin brightnesses are plotted as a function 
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of projected distance from the nucleus. The brightnesses are shovm for the Pi 
transition of the IR (1,0) band using two different excitation models: the 
circles, 0, are the non-equilibrium (DGEAR) results and the triangles, A, 
represent the steady-state solution. Significant differences between the two 

3 

models occur only when the projected distance from the nucleus is between 10 

4 

and 10 km. The spatial distribution of the brightness in the inner ccma is 

different from that of the column density, which varies inversely with the 

2 8 

projected distance fron the nucleus. The CO production rate is 2x10 mol 

-1 29 -1 

s , the H^O production rate is 2x10 mol s , and the heliocentric distance 
is 1 AU, 

Figure 3. Colimn densities and colunn brightnesses are plotted as a function 
of projected distance from the nucleus. The brightnesses are shown for the PI 
and P7 transitions of the IR (1,0) band. The colimn density varies inversely 
with the projected distance fron the nucleus out to v'10^ km. The PI and P7 
brightnesses have different spatial distributions, with neither being directly 
proportional to the column density in the inner coma. The col twin brightnesses 

have been calculated using the non-equilibrium (DGEAR) model. The CO 

28 “1 29 "1 
production rate is 2x10 mol s” , the H 2 O production rate is 2x10 mol s” , 

and the heliocentric distance is 1 AU, 

- 1-1 '' 

Figure 9. Average g-factors (in photons s mol ) in a 1 ' field of view for 
the IR (1,0) band when the geocentric distance is 0.62 AU. Results from three 
different excitation models are shown: the circles, 0, represent the 

noivequilibri jm (DGEAR) solution, the triangles, A, are the steady-state 
results, and the crosses, +, represent the distribution when rotational LTE is 
achieved in the upper vibrational level (v=1). Notice that the intensity in 
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the IR (1,0) band Is distributed anong many lines in this example. The CO 
production rate is 2x10^® mol the H^O production rate is 2x10^^ mol s”\ 

and the heliocentric distance is 1 AU. 

A 

Figire 10. Same as Figure 9, except for a 4” field of view. In this case the 
intensity in the band is concentrated into the low J transitions. The largest 

A 

g- factor here is almost three times the largest value attained for the 1’ 
field of view. 
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II. RESULTS CF THE CALCULATION FOR 
STANDARD BRIGHT COMET, 

CCWET IRAS-ARAKI-ALCOCK, 

AND COMET HALLEY 
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INTRODUCTION 


The results of the non-equilibrium calculation for different comet 
parameters are presented in this section. The details cT the calculation are 
given in Part I. Here, we give results for three different comets; a standard 
"bright” comet, comet IRAS-Araki-Alcock (1 983d), and comet Halley. 


28 —1 

For the bright comet, we use a CO production rate, of 2x10 mol s" , 


a water production rate, ^ solar flux at R=1 AU. 

To gauge the effects of solar pumping on the CO excitation, we vary the solar 
flux to approximate the conditions of the comet at R=2 AU and R=0 " AU 
(production rates and other parameters remain Identical to their val . R = 1 
AU). For comparison, we also use the standard bright comet parameters to 
calculate a steady-state solution and a non-equilibrium solution with an 
isothermal temperature profile (T=200K) throughout the coma. 


For comet IRAS-Araki-Alcock (1983d), we use H2O production rates and upper 


limits on the CO production rate from lUE measurements (Feldman, private 
communication). We calculated two cases for comet 1 983d; IRAS1 with 


,27 


-1 


,26 


Q^q=1x 10 mol s and <3^20 IRAS2 with Q^q=5x 10 mol 

-1 


5 and Qu„=40xQ , and R = 1 AU. 
n20 CO 

Three different observing aspects for comet Halley are used to calculate 


model predictions; HALLEY1 with Qj,q=4.3x 10^'^ mol s”\ Sl20" ^ = 


p O 1 

and geocentric distance (a) of 0.62 AU, HALLEY2 with Q_ =1.3x10 mol s” , 

C 0 


‘^20'"'°’"^C0’ HALLEYS with Q^^=5.7 x1 0^*^ mol s"\ 


20 -“' ^=1.33 AU, a=0.42 AU. 


‘CO 


Ki 


Results for the standard bright comet, comet IRAS-Araki-Alcock, and comet 
Halley are given in four groups: 

Group I shows contour plots of the fractional populations for CO 
rotational transitions and excitation tenperatures as a function of distance 
from the comet nucleus. 

Group II shows the fractional population and excitation temperatures of 
selected rotational states of CO as a fuction of distance from the nucleus. 

Group III shows coltmn densities and brightnesses of selected CO IR (1,0) 
tr an sit ions . 

Group IV. shows aperture averaged intensities and g-factors for the CO IR 
(1,0) band using a and field of viaw. 


Figure Captions 


NOTE: Since there are flgires, captions cannot be given for each. The 

figures have been divided Into four groups. A general figure captions Is givsn 
for each group followed by a listing, with a short id«itlflcatlon for all the 
figures in that group. 

I. Contour plots showing relative populations of CO rotational states and 
excitation temperatures for CO rotational transitions as a function of distance 
from the nucleus. The plots are arranged in pairs consisting of the plot for 
the relative populations followed by the excitation temperature plots. 
Different pairs of plots are distinguished by different comet conditions and/or 
the use of different excitation models. For the relative population plots, each 
vertical line represents r different rotational state, starting with J=0 at the 
extreme left and incrementing by one until J=20 is reached at the extreme right. 
The contours are given in logarithmic units so that, for exanple, the cirve 
labeled -1.0 represents a relative population of 10%, Excitation temperature 
are calculated between adjacent rotational levels (see text for explicit 
formula) with each transition represented by a vertical line. Excitation 
temperatures for the lowest frequency transition (J=1-»0) are given along the 
left most dotted vertical line. Regions in which a transition becomes inverted 
(negative excitation temperatures) are cross-hatched. In all cases, plots are 
shown for both the v=0 and v=1 vibrational levels of the ground electronic 
state. Also, the regions between the vertical lines have no physical 
significance. Individual figure identifications follow. 
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lal 

and 

Ia2: 

Ibl 

and 

Ib2: 

Id 

and 

Ic2: 

Idl 

and 

M2: 

lei 

and 

Ie2: 

If 1 

and 

If2: 

Igl 

and 

Ig2: 

Ihl 

and 

Ih2: 

111 

and 

Ii2: 


and 

Ij2; 


28 

DGEAR (nory-equilibrium) model with Q^qs2x 10 mol s” *^20"^^*^C0' 
R=1 AU. 

Steady-state model vdth Q^q=2x 10^® mol s"\ Si20"^®’^CCj' 

DGEAR moda with Q^qs2x 10^® mol s“\ Sl20"^®*®C0* 

28 --1 

DCEAR model with Qj.p=2xl0 mol s , AU. 

Sane as (a) except using an isotherm^ cana with T, . s200K. 

kinetic 

27 -1 

DGEAR model for the conditions of IRAS!: Q^^slxlO mol s , 

*''• 

Pfi —1 

DGEAR model for the conditions of DIAS2: Q^^=5x!0 mol s” , 

Sl20=''‘’'‘«C0’ 

DGEAR model for the conditions of HALLEY1: Q _=A,3xlO^^ mol s”\ 

DGEAR model for the conditions of HALLEY2: Q^^=1.3x10^® mol s*\ 
S{20''^°’'^CC’ «=0-88 AU. 

DGEAR model for the conditions of HALLEY3: Q_.=5.7x10-^ mol s"\ 

w V 

O^2o’^®*^r0' ^* = ‘'•33 AU. 


II. The fractional populations of selected rotational states of CO and the 
excitation temperatures of selected rotation^ transitions are plotted as a 
function of distance from the nucleus. The figures are divided into five groups 
of twelve plots each. Each group represents a different set of comet conditions 
and/or the use of different excitation models. The first half of each group 
(six plots) shows fractional populations for the J =0, and 10 rotational 
states, for both the v=0 and v=1 levels. In these figures, the solid curve 
represents the fractional population achieved assuming rotational LTE. The 
second half of each group shows excitation temperatures for the Js1-K), 5^4, and 
10-^ rotational transitions, for both the v=0 and vsl levels. For these plots. 
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the solid 
inversions 
Ir^ iv id ual 

II a1-a6: 

II a7-12: 

II b1-b6: 

II b7-b12; 

II c1-c6; 

II c7-cl2: 
II d1-d6: 

II d7-d12: 
II e1-e6: 


II e7-e12; 


curve is the assined kinetic temperature profile and population 
(negative excitation temperatures) are indicated by arrows, 
figure identifications follow. 


DGEAR (norwequilibriUB, and Steady-State fractional populations. 
Q^q=2x 10^® mol s"\ 

26 

DGEAR and Steady-State excitation temperatures. Q^^^xlO mol 

s”\ Si20”'°*®C0* 

Fractional populations using DGEAR model for conditions of: 

IRASl: Q^q=1x 10^ mol s"\ 

IRAS2: 

Excitation temperatures using DGEAR model for conditions of IRASl 
and IRAS2. 


Fractional populations using DGEAR model for R=2 AU and R=0.5 AU. 

28 —1 

In both cases, Q^j^=2xl0 mol s and 

Fractional temperatures using DGEAR model for R=2 AU and RssO.S AU. 

Fractional populations using DGEAR model for conditions of; 

HALLEY1: Q^q=4.3x 10^^ mol s"'* , Q„20^’°’^C0 

28 1 

HALLEY2: Q^,q = 1.3x 10 ° mol s" , 

Excitation temperatures using DGEAR model for conditions of 
HALLEY 1 and HALLEY2. 


Fractional populations using DGEAP model but with two different 

kinetic temperature profiles. One case employs the Marconi and 

Mend is (1982) kinetic temperature profile (same as used above in 

all cases), while the other assitnes an isothermal cona with 

T. , . , =200K. 

kinetic 

Same as II e1-e6, except for excitation temperatures. 
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III. Column densities and column brightnesses are plotted as a function of 
projected distance frcm the nucleus. 

Ill a-c: These figires compare results for the DGEAR (non-equilibrium) and 

Steady-State models for the PI, P6, and P8 lines, respectively, of 

the CO IR (1,0) band. For the PI transition, significant 

differences between the two models occur only when the projected 

3 4 

distance f»“om the nucleus is between 10 and 10 km. For the P6 
and P8 lines differences can also occur for small projected 
distances (<10^ km), but for these lines the differences between 
the two model results are generally small (<101). 

III d-f: These figures are the results from the DGEAR model comparing 

spatial brightness distributions for two different lines. Figure 
Illd compares PI and P2, Ille compares PI and P4, while Illf 
conpares PI and P7. For all of the lines in the CO (1,0) band, 
the spatial brightness distribution in the inner coma is different 
fron that of the column density, which varies inversely with the 
projected distance fhora the nucleus. 

28 “1 

All of the above results have been calculated using Q^q=2x 10 mol s” , 

IV. Aperture- aver aged intensities and g-factors for the CO IR (1,0) band are 

A A 

shown for two fields of view (4" and 1'). Individual figure identifications are 
given below. 
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IVa: 

IVb: 


IVc: 


IVd: 
IVe: 
IVf : 
IVg: 


DGEAR ( ''on-equilibriun) and Steady^State model results using 
Q^q=2x 10^® mol s”\ Si20’^°*®CO’ 4=0.62 AU, and a field of 

A. 

View (FOV) of 4". 

A 

DGEAR model using FOV^" for the conditions of: 

HALLEY1: Q^q=4.3x 10^ mol s*'\ Si2o“^°*W ** = ^-53 AU, 4=0.62 AU 

pQ ^1 

HAaEY2: Q^q=1.3x10‘^ mol s , ^ao”^°^CO' 4=1.0 AU 

HALLEY3: Qj.q=5.7 x10^^ mol s‘\ ^ = ^*33 AU, 4=0.42 AU 

A 

DGEAR model using FOV =4" for the conditions of: 


BRIGHT: 

Q^O=2.0x10"» 

mol 

-1 

s , 


R = 1 AU, A=0.062 AU 

HAaEY2: 


mcl 

-1 

S , 

(^2o=’0*Qco* 

R=0.88 AU, 4=1.0 AU 

IRAS1: 

Qco=l.0x10=^ 

mol 

-1 

S , 

Sl20''°**^C0’ 

R = 1 AU, 4=0.041 AU 


Same as (a), except using F0V=1 ' 

A 

Same as (b), except using F0V=1' 

A 

Same as ( c) , except using F0V=1 ’ 


g- factor X usir^ FOV =4", 4=0.62 AU, 

and 

the 

DGEAR 

model 

for 

three different values of R: 1 AU, 0,5 

AU, 

and 

2 AU. 

In 

all 

28 —1 

cases, Q^q=2x 10 mol s“ and 

A 

Same as (g), except using F0V=1'. 

♦ 
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Figure Illf. 
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Figure IVc. 
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Figure IVg. 
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